Effective theory of excitations in a Feshbach resonant superfluid 
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A strongly interacting Fermi gas, such as that of cold atoms operative near a Feshbach resonance, 
is difficult to study by perturbative many-body theory to go beyond mean field approximation. Here 
I develop an effective field theory for the resonant superfluid based on broken symmetry. The theory 
retains both fermionic quasiparticles and superfluid phonons, the interaction between them being 
derived non-perturbatively. The theory converges and can be improved order by order, in a manner 
governed by a low energy expansion rather than by coupling constant. I apply the effective theory 
to calculate the specific heat and propose a mechanism of understanding the empirical power law of 
energy versus temperature recently measured in a heat capacity experiment. 



In the last few years, the rapid development of ex- 
periments of atomic Fermi gas using Feshbach resonance 
has revitalized theoretical interest of strongly interacting 
many-body system. Recent experiments, for instance, 
the collective excitation J„, and heat capacity point 
to the need of a dynamical theory beyond mean-field 
BEC-BCS crossover theory. A basic difficulty with the 
crossover theory is that perturbative expansion is not 
justifiable if going to next order for strong interaction. 
Various efforts have been made on the basis of hydro- 
dynamic theory 0, 0, IS 1 which in general explain 
the frequency of collective modes particularly well on the 
resonance point. A disadvantage of hydrodynamic the- 
ory is that fermionic excitation is hard to incorporate 
beyond the density and velocity field approximation. In 
this paper, I develop an effective field theory for the res- 
onant superfluid that retains both fermionic and bosonic 
excitations explicitly based on broken symmetry. An ad- 
vantage of this theory over a pure bosonic theory is that 
it allows one to study, for example, the damping of col- 
lective modes due to scattering with fermionic quasipar- 
ticles. The present approach is thus different from and 
complimentary to the existing crossover and hydrody- 
namic theories. 

Effective field theory of excitations. The physical sys- 
tem we are interested in is a ultra-cold gas of fermionic 
atoms (such as ^Li or ''"K) with two internal (hyperfine- 
spin) states in a broad Feshbach resonance. It is known 
from both theory and experiment (see, for instance, 9] 
and ^3 J respectively) that the molecule fraction is small 
from the resonance regime through the entire fermionic 
side. The microscopic theory is then described by the La- 
grangian (in the path-integral framework with imaginary 
time, T — it): 



L 



where ■00- {x) describes fermionic atom fields at a: = (x, t) 
and spin state a =t, i, is the coupling constant defined 
by 5 = (a the s-wave scattering length) , and ^ is the 
chemical potential. (Throughout, summation is implicit 
over repeated spin indices; 7i = 1 is set in units.) 



For attractive interaction g, one can perform an exact 
Hubbard-Stratonovich transformation to introduce the 
auxiliary pair fields, A(x), and its complex conjugate. 
Then, the Lagrangian becomes 

+(V'|VIA + c.c.) + i|A(x)p}. (2) 

Up to this point, no approximation has been made. One 
may proceed then to derive the gap and chemical equa- 
tions, following the BCS-BEC crossover theory procedure 
(for instance, Ref. ^3 ^'^d references therein). But we 
are interested in the dynamics of low energy excitations 
about the superfluid state. 

The atomic gas is in strong coupling regime, having 
a scattering length greater than or comparable to the 
Fermi wavelength kp"^ {\kFo\ ^ !)■ Conventional per- 
turbation theory does not directly apply, at least not in 
a standard, controllable manner. To have a controlled 
theory for the description of quantum dynamics beyond 
mean field theory, one may seek for available different 
approach. 

The broken symmetry of the superfluid state provides 
such an approach — effective field theory ■ The atomic 
gas conserves the total number of particles. This conser- 
vation law corresponds to a U(l) symmetry of the phase 
transformation. For the above model, the U(l) symme- 
try transformation is 



'A(x), 



(3) 



with a an arbitrary, constant phase. 

Recall that the superfluid state is a Bose-Einstein 
condensate of fermionic atom pairs, which breaks the 
U(l) symmetry. The order parameter acquires a finite 
value below the superfluid phase transition temperature, 
(A(x)) = Aq 0. Fluctuations about this superfluid 
ground state are excitations. 

Here I derive an effective theory of the phase fluctua- 
tion, directly based on the broken symmetry. We start 
by treating Aq — the groundstate or thermal mean value 
of the order parameter — as an input parameter. It may 
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be taken directly from experiments; of course it can also 
be taken from the result of Monte Carlo, or any mean 
field crossover theory but that will require further jus- 
tification. We are concerned with phase fluctuations in 
space and time, so write the order parameter field 



(4) 



(The factor 2 is inserted to make it explicit that the pair 
field carries two units of mass.) We have set to ignore 
the amplitude fluctuation for the purpose of simplicity; 
this is justifiable for it is known gapped and, as can be 
checked below, does not directly coupled to phase except 
via fermions. Instead of making a Gaussian expansion — 
a usual practice, we transform non-perturbatively the 
fermion fields at each space time point as follows: 

Mx) = M^y"^''^ , i'Ux) ^ rA^)e-'"^^K (5) 

The transformation is designed to eliminate the phase 
fluctuation 0{x) dependence from the off-diagonal pairing 
potential terms, {ip*tfj*A + c.c). As a result, the 9{x)- 
dependence can only arise from the kinematic terms of 
the fermion theory sector. Collecting all terms of tp and 
9, we have 



§^)7^, -f (AoT^xV^T + c.c.) (6) 



where D^, = + iduO{x) , ly = T,x,y, z . The U(l) sym- 
metry guarantees the phase field to appear only through 
derivatives, since energy cannot depend on a constant 
phase. The form of coupling between the phase and 
fermion fields is exact from the model. This theory is no 
perturbation, apart from excluding the amplitude fluctu- 
ation. 

For those who are familiar with gauge theory, the time 
and spatial derivatives of 6{x) together behave like a four 
component gauge field. After the transformation , the 
fermions ^pa see a constant off-diagonal pairing potential 
A{x) = Ao but interact with the (Goldstone) phase field 
0{x) by a means similar to 'gauge' coupling. 

Landau superfluid hydrodynamics. The effective the- 
ory jni has a very strong implication for superfiuid hy- 
drodynamics. To reveal that, let us collect all terms of 
first order in 9 from Jif^g'. 



+fi{x)idr9{x) +W0{x) ■ 3{x) 



(7) 



where n{x) = ^|J*{x)^pa■ix) = ip*{x)'ijja-{x) and 3{x) = 
— 2^(?/'JVV'ct ~ V'0i'0CT) are the fermion density and cur- 
rent operators. The first term implies that the fermion 
total density and the order parameter phase are canoni- 
cally conjugate. The amplitude fluctuation of the order 
parameter is separate. [This can be checked by shifting 
Ao by a quantum fluctuation A'(x), Aq ^ Ao + A'(x).] 
To derive the hydrodynamic theory, we add to the 
term 



A^ 



(8) 



and do path integrals over the new flelds p{x), X{x). This 
is allowed because the path-integral over X{x) gives a 6- 
functional, which sets the equation 



p{x) = ^Ij*{x)'iI;„{x) , (the density held) 



(9) 



and the path integral over p{x) gives a constant. The 
equation justifies p{x) as fermion density field. The 
added term is designed to extract the effective den- 
sity field. Up to now the transformation of is exact. 
Performing path integrals over fermion fields produces 
the effective action: 



S — SpB + S\ + 



(10) 
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+ J^^, ^^jix - x')d,Mx)d^,eix') 

Sx = I ip\ - \ 9(x - x')\(x)\{x') (11) 

where '. . . ' stands for higher powers of 9 and/or A fields. 
Here, '3 and are the density and current correlation 
functions defined in the homogeneous ground state (so 
the subscript '0'): 

^{x-x')^~{h{x)n{x'))a and (12) 
S^^{x-x') = -i^X,,[x)Jj{x'))Q, i,j^x,y,z. (13) 

This is a low energy (derivative) expansion of the collec- 
tive fields, rather than a perturbation in some coupling 
constant. The derivative nature of X{x) is less obvious, 
but immediately becomes manifest when examining the 
quantum equation of motion. 



6p 



Therefore the power expansion on A is equivalent to a low 
energy expansion. The field A is auxiliary in nature. We 
can calculate the effective couplings of 9 and p mediated 
by A(a;) perturbatively in low energy limit. The result- 
ing effective action is the equivalent of Landau superfluid 
hydrodynamics: 



^'pfi — 



So 



x') ^p{x)p{x') 



(14) 



where Spe is the same as given above. 

Effective theory of phase fluctuations. Single fermion 
excitations are gapped at all energies below Ao. The 
Goldstone boson 9{x) interacts with fermion through 
particle-hole pairs. Therefore, for collective excitations 
of energy below 2Ao, we can integrate out all fermionic 
degrees of freedom and have an effective theory for 9 
alone. Starting from the theory ((HJ, we write it into a 
quadrature of fermion fields in the Nambu spinor space, 



d^xdr 



("0* (a;)>i(a;) 



fi^^ix'r 



(15) 
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with 

i^x.x' = [if"(x)+i^''(a;)](53(x-x')<5(T-T'), (16) 
K'^ix) = Wr - T3(|^ + /i) + TiAn , (17) 

K'{x) = +r,^^^e, (18) 

where {1, ri,r2,T3} are the identity and Pauh matrices 
for the Nambu spinor space. Integrating out fermion 
fields gives the effective action for the phase mode Sg = 

After some straightforward calculation, the effective 
action of the theory is found to be (in momentum- 
frequency space) 



Sf) = 



L0' + 



^l|('?) „2 , «o.„2 
2 4 ^ 2m4 



0{qrO{q) 



1 

+o( 



>(2A^' (2A^) (l^igher powers of 9) (19) 
where ng is the total fermion density and 



}- ^ tr{T3G(p)T3G(p + q)} , q = (q, it^) 



/3y 



('^) = ^ E Pitr{G(p - f )G(p + |)} , II q. 



G{p) 



with er 



Tsep + riAo 



[iuj' - Ep){iuj' + Ep) 
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fermionic quasiparticle spectrum. We are working in the 
imaginary time formalism: uj and oj' are the Matsubara 
frequency for bosons and fermions, respectively; f3V= 
space-time volume; and G{p) is the fermion propagator. 
For readers who are familiar with diagrammatic calcula- 
tion, &{q) and ^{q) are directly related to the density 
and current correlation functions defined in Eqs. p2(l and 
(|13|l by Fourier transformation. 

One can directly evaluate &{q) and S^{q)- The cal- 
culation is tedious but otherwise elementary. At zero 
temperature T — Q, the result is 



Ag. Bp is the 



&{q) 
^11 (g) 



1 



1 



{iuj -E+- E^){iLU + E+ + E^) 



,(20) 



where J^ = J ^^^^3 , and the subscript ± stands for e± = 
£p±q/2 and likewise for E±. (Note the important sign 
difference in the 'coherence factor' between ^ and 3^.) 

Here to an important point regarding the perturbative 
expansion. The expansion is sorted by the powers of mo- 
mentum and frequency of the field ^(q, iw). Unlike the 
Gaussian fluctuation theories, the present approach does 
not require the phase fluctuation field itself small but 
slowly varying in space and time. This approach is also 



quantitative: The low energy and momentum scales are 
set by Aq and Aq/vf, respectively; One can improve the 
theory order by order and estimate the next order cor- 
rection. A more complete development of this approach 
will be given in a future study. The essential point of why 
this approach works is that the phase fluctuation (Gold- 
stone) field must appear through derivatives dictated by 
the broken U(l) symmetry. 

For those familiar with field theory, I should be careful 
about approximation in the expansion of H19|l . The gen- 
eral form is no approximation, but the coefficients & and 
however, could not be derived in any 'honest' calcu- 
lation for strongly interacting gas if we had included the 
amplitude fluctuation too. So the coefficients should be 
treated as good as mean field approximation. Nowadays, 
a more typical application of Weinberg effective field the- 
ory in particle physics is to directly construct ingredients 
of proper symmetry transformation, write a general form 
of effective Lagrangian, and match the unknown coeffi- 
cients with experiments rather than derive. We will not 
discuss that general approach here. 

The ujq <C 2Ao limit. This corresponds to the limit 
of long wavelength, q ^ 0. Expanding the and ^ of 
Eq. H20|) in powers of q^, we get the leading order 



S){q = 0) = - 
^{q = 0) = 



d'p Al 
(2^)3 i?3 



A'^mkp 



J (21) 
(22) 



where A = Ao/ef, M = m/cf and -JP = 
f„°° dx- — ^ .„,„ . [ef is defined as the Fermi energy 

at the non- interacting limit.] The difference between ^ 
and can be seen from the coherence factor in their 
expressions. In the case of S^{q), the factor is exactly 
vanishing at q = 0. 

Then the effective action of the phase mode becomes 



^ 



no „2 
m. 1 



(2Ao)'" (2Ao)* 



(23) 



The energy spectrum of the collective mode is obtained 



by analytical continuation iuj 



iO^ and by finding 



the zero of the coefficient of the quadratic 0-term. The 
spectrum is 



Vs\q.\+0 



(24) 



where Vs is the sound velocity. 

For small Aq compared with the Fermi energy ep (BCS 
limit), y « ^[1 + 0{Al)]. Using the relation kp = 

(Btt^/io)^/^ and vp — kp/m, this reproduces the sound 
velocity Vg — vp/\/3 — the well-known result in the BCS 
limit obtained in the theories of Bogoliubov, Anderson, 
and others later on [lllll3| . 
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Specific Heat. Recently, Kinast et al |3] reported a 
measurement of the heat capacity of a strongly inter- 
acting Fermi gas of ^Li atoms, slightly detuned to the 
BCS side from the resonance. The heat capacity reveals 
a kink at a temperature about T = 0.27Tp, arguably 
an evidence of superfluid phase transition. While the 
crossover theory employed in the paper agrees with the 
data in general, it deviates from the data noticeably at 
low temperature below the transition, as shown in their 
figure of energy input versus temperature. 

As an example of application, I will use the effective 
field theory (|23|l to calculate the specific heat for tem- 
perature well below the superfluid critical temperature 
Tc {T <^ Tc ^ 2Ao/fc_B). All fermionic quasiparticle 
excitations are gapped, so the only contribution to the 
entropy must be due to the collective superfluid phonons 
with a linear dispersion. The thermal energy per unit 

volume due to the phonons is ^ = ^ 30^^'' • ^he coef- 
ficient factor (including Vg) is model dependent; but the 
power law is universal in the limit of T <C Tc. Dif- 
ferentiating the thermal energy gives the specific heat, 
c = — — . By measuring the coefficients of the T"^ 
law, one can determine the sound velocity and check the 
theoretical result. Fig. ^ shows how Vg and c vary with 
Aq between weak and strong coupling limits — a result 
seeming not reported anywhere else yet. 




FIG. 1: Sound velocity and specific heat in a reso- 
nant fermionic superfluid. The chemical potential is self- 
consistently determined when varying Aq with density fixed. 
Parameters: kp — (Sn^no)^^^; ep = kp/(2m); vp — kp/m. 
The BEG limit (molecular side) is plotted to just show the 
trend; otherwise it is beyond the scope of validity of the 
present theory. 

The power law of the energy differs from the em- 
pirical E cx T^-^^ 5]. Bulgac [lj| examined the role of 
collective excitations by focusing on the unitary regime 
and predicted a complicated T-dependence in different 
temperature regime, from exponential to i? cx T^. Here 
I provide an explanation other than those discussed in 
Ref. H Q. Away from the limit of T < T^, we 
now move to the temperature regime slightly below Tc 
(0 <C T ^ Tc). There can be preformed Cooper pairs but 
the fraction of condensed Cooper pairs is small. In our 



effective field theory description, that corresponds to a 
small amplitude of order parameter, Ag. The superfiuid 
phonon still exists but looses the linear dispersion char- 
acteristic. This fact is actually implicit in Eq. H24I) for 
which the expansion ought to include nonlinear powers 
of |q|/Ao when Aq is relatively small. In fact, for tem- 
perature near Tc, the expansion should be done another 
way around, in terms of small Ao/|q| instead. At the 
end, the superfluid phonon and non-condensed Cooper 
pair, coupled together, all are expected to have a disper- 
sion quadratic in momentum, Wq ~ q^. Such a Bose 
gas gives a temperature dependence as i? cx T^/^ by 
simple power counting. For strongly interacting gas, 
fermionic excitations are believed to have a pseudogap 
far greater than Aq — sometimes called coherent super- 
fluid gap. Therefore fermionic excitations do not alter 
the above power law. (There might be contributions due 
to gapless fermions on the trap edge if the effect of trap 
were included, but that seem not as important in the 
thermodynamic limit.) In summary, the power law is 
for T/Tc -> and is argued to be T^/^ for < T < Tc. 
Perhaps, in my opinion, an interpolation of these two 
limits can flt the data better, entirely consistent with the 
empirical T'^-^^ law from a single-exponent fit. 

All results apply on the fermionic side of a broad Fesh- 
bach resonance, when Aq > (the larger the better, con- 
trary to weak perturbation theory) . One may easily gen- 
eralize our approach to develope a more complete effec- 
tive theory valid for the molecular side as well, by starting 
with a Bose-Fermi resonance model of both atoms and 
molecules. 

Finally, I learned that the low energy expansion was 
used long ago to derive time-dependent Ginzburg-Landau 
theory for a BCS superconductor by Abrahams and 
Tsuneto |l5l |. But their expansion was done in a very 
different manner and assumed the (weak coupling) BCS 
limit which was quite alright at that time. Another use- 
ful difference is that the present effective theory, having 
retained fermions explicitly, is convenient to study, for 
example, the damping of collective modes due to decay 
into (fermionic) quasiparticles. 

I am grateful for discussions with D. Son, J. Thomas, 
and F. Zhou. I thank the Institute for Nuclear Theory at 
the University of Washington for its hospitality through 
the Program on Quantum Liquids and Gases and the De- 
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pletion of this work. 
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